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Abstract 

We consider the logistic map over quaternions H ~ R'* and different 2D projections of 
Mandelbrot set in 4D quaternionic space. The approximations (for finite number of iterations) 
of these 2D projections are fractal circles. We show that a point process defined by radiuses 
Rj of those fractal circles exhibits pure 1// noise. 
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1 Introduction 

1// noise is observed in large diversity of real life and artificial systems, which behavior is usually 
defined by a complex interaction of many components. Complexity of the system usually assumes 
that long-term correlations are observed. Examples are processes and experimental data in con- 
densed matter, traffic flow, quasar emissions, music, biological and medical systems, economic and 
financial data, human cognition and even distribution of prime numbers (see PP and references 
herein). 

Fluctuations of signals defined by time series obtained from such systems are found to be cha- 
racterized by a power spectral density S{f) diverging at low frequencies / like 1//", here a is some 
real parameter. 1/f (a « 1) noise is an intermediate between the white noise {a — 0) with no 
correlation in time and the random walk (Brownian motion) noise (a — 2) with no correlation 
between increments. Note that Brownian motion can be obtained integrating white noise and 
that taking the integral of the signal increases the exponent a by 2 while the inverse operation of 
differentiation decreases it by 2. 

Parameter a is closely related to the Hurst exponent H . It is known that fluctuations which 
are fractionally homogeneous, i.e. unifractal or uniscaling, can be quantified by a single coefficient 
H and a single exponent a 

Possible generalization leads to multiscaling or multifractals, with the exponent H dependant 
on time. Therefore multifractal processes are characterized by a set of scaling relations or power 
laws with correspondingly many exponents a [21 . 



2 Point processes and 1// noise 

In many cases, the intensity of some current can be represented by a sequence of random (however, 
as a rule, mutually correlated) or pseudo-periodic pulses A-k[t — tk). Here the function Ak{(p) 
represents the shape of the fc-th pulse having an inffuence to the current I{t) in the region of transit 
time tk. The intensity of the current in some space cross-section may, therefore, be expressed as 

i{t) = J2Mt-tk). 

k 
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It is easy to show that the shapes of the pulses mainly influence the high frequency, / > Atp 
with Atp being the characteristic pulse length, power spectral density of I(t) while fluctuations of 
the pulse amplitudes result, as a rule, in the white or Lorentzian but not 1// noise. 

Therefore, we restrict our analysis to the fluctuations due to the correlations between the transit 
times tk and hence we can replace the function Ak{t — tk) by the Dirac delta function S{t — tk). 
The current (see Fig.P) is then expressed as 

m^Y.^{t-tu). (1) 

k 

Following this approach, instead of current I{t), we further deal with point process^ defined by 
the sequence ti,t2, ■ ■ ■ ,tjq , . . .. 



I(t) 
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Figure 1: Current I{t) vs time t defined by formula. Such dependencies appear when registering 
the consecutives heart beats, cars on a highway passing through the reference point, transactions 
in financial markets etc. 

The power spectral density of the current is defined as 

2 

(2) 

where ijv] is assumed to be the interval of observation. 

In this approach the power spectral density of the signal depends on the statistics and cor- 
relations of point process (the transit times tk) only. It is well known that sequence of random, 
Poisson, transit times generates white (shot) noise, for example. 

In ,4^ we proposed simple analytically solvable model for producing point process resulting in 
S{f) ~ 1// (a = 1) noise. Discussion on the origin and universality of 1// noise was continued 
in [SIE]. Some further work, related to the applications of the theory of point processes and 1// 
noise to econophysics, was done in |71|S]. 

3 Quaternions and other hypercomplex numbers 

Complex numbers C ~ R^, along to their real predecessors R, are widely used in nowadays math- 
ematical modeling and scientific computing. Beside others, they have important applications in 
theories of complex systems, fractals and signal processing: famous Mandelbrot and Julia fractal 
sets are defined in C, spectrum (Fourier transform) is defined as integral of complex function etc. 

There are some clues that we should not stop with the computations in R and C, and that 
further generalization to quaternions H ~ R"* (introduced by Hamilton) or even octonions G ~ 
(introduced by Graves) are particulary interesting and valuable, even though the role of these 
hypercomplex numbers is not widely understood yet. 

In order to define hypercomplex algebras, one has to consider not only two algebraic operations 
-t- and X , but also one geometric map: a; i— > x, where x denotes the conjugate vector of x. 

The three operations are defined recursively as we define the algebras, in the following manner. 
Let Ak be the real hypercomplex algebra of dimension 2^ , k >1. It is constructed recursively as 



S{f) ^ hm - 
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k=l 
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Ak = Ak-i X Ak-i by means of the three following operations: 



addition: (a, b) + (c, d) 

conjugacy: (a, b) 

multiplication: (a, b) x (c, d) 



{a + c,b + d), 
{a, -b), 

(ac — db, da + be) , 



where ac denotes a x c in A^-i- For k = 0, Ao is taken to be the field M with the arithmetic 
operations + and x , the conjugacy map being the identity on R: a ^ a = a G M. This construction 
is known to algebraists as the Cayley-Dickson doubling process. 

About computations with hypercomplex numbers, and why only real numbers, complex numbers, 
quaternions and actions are suitable for computations see (S) ll(J| and references herein. 

Explicitly multiplication in H can be expressed as {a,b,c,d) x {a' ,b' , c' , d') = {a" ,b" , c" ,d"), 
with 



b" = ab' 



c 
d!' 



ac 
ad! 



ba! - 
ca' - 
da' 



cd' - dc' 
db' - bd' 
- be' - cb' 



4 1// noise in quaternionic Mandelbrot set 

We consider the logistic map over quaternions H ^ R** 



0,1,. 



(3) 



with given initial value zo, for example zq = (0.5, 0, 0, 0). The logistic map Q has been extensively 
studied over R (real numbers) and C (complex numbers). Despite its great simplicity this map 
exhibits an extremely complex behaviour. The study of (j^J on R gives birth to the Feigenbaum 
tree while the analysis of (j3Jl on C leads to the famous Mandelbrot and Julia fractal sets. 

Further we deal with 2D projections of Mandelbrot set in 4D quaternionic space. Any two 
components of r are set to zero, while the remaining two vary. For example, 

Mi2 = \ {ri,r2) ■■ r = (ri,r2,0,0), lim \zk\ < oo L 



Mo 



(r2,r4) : r = (0,r2,0, r4), lim \zk\ < oo 

k — >oo 



Note that Aii2 is just the famous Mandelbrot set in 
and M23 = M2A = ^^34- 



We also show that M12 — M13 — M14 




Figure 2: Approximation (after 50 iterations) of Mandelbrot set A^23 (one gets exactly the same 
for M24 or M^i). 

The approximations (for finite number of iterations) of Mandelbrot set A^23 = A^24 = ^^34 
(near its boundary) are fractal circles (see Fig.O, dependant only on radius R ~ \/^\ + ''3 ■ 
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Figure 3: The number of iterations needed to reach \zk\ > 10^" plotted vs radius R when computing 

M23- 
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Figure 4: Power spectral density S{f), defined by l@J, vs frequency / with N = 796474. The plot 
is compared to the function 1//. 

Define the point process Rj as the values of radius of each circle - mathematically they are the 
values of R, small change of which result in significant change of number of iterations needed for 
\zk\ to reach "infinity" (10^° for example). The values Rj correspond to peaks in Fig.|2| 

According to ^ , the power spectral density of such point process is defined as 



S{f) 



Rm — i?1 



N 



(4) 



here N is the volume of point process data {N — > cx) as Rj recording resolution increases). 

We obtain (see Fig.0J that S{f) ^ 1//, i. e. radiuses Rj of fractal circles in Mandelbrot set 
A^23 exhibit pure 1// noise (a = 0) or unifractal noise. 
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